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4^ . ABSTRACT 

Oh- 

Context. The Parkes Galactic All-Sky Survey (GASS) is a survey of Galactic atomic hydrogen (H i) emission in the southern sky 
observed with the Parkes 64-m Radio Telescope. The first data release was published by McClure-Griffiths et al. (2009). 
Aims. We remove instrumental eff'ects that aff'ect the GASS and present the second data release. 
^ ■ Methods. We calculate the stray-radiation by convolving the all-sky response of the Parkes antenna with the brightness temperature 

' distribution from the Leiden/Argentine/Bonn (LAB) all sky 21 -cm line survey, with major contributions from the 30-m dish of the 

Instituto Argentino de Radioastronomia (lAR) in the southern sky. Remaining instrumental baselines are corrected using the LAB 
^vq ■ data for a first guess of emission-free baseline regions. Radio frequency interference is removed by median filtering. 

^ ■ Results. After applying these corrections to the GASS we find an excellent agreement with the Leiden/Argentine/Bonn (LAB) survey. 

\^ ' The GASS is the highest spatial resolution, most sensitive, and is currently the most accurate H i survey of the Galactic H i emission 

QQ \ in the southern sky. We provide a web interface for generation and download of FITS cubes. 
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1 . Introduction ISM in the soft X-ray energy range below 1 keV. The attenuation 

of the X-rays by the Galact ic interstel lar medium increases with 

Atomic hydrogen (Hi) is ubiquitous, it traces the interstellar decreasing X-ray energy (^|2TO3). The emission of active 

medium (ISM) m disk galaxies over a broad range of physical Galactic nuclei at high redshifts (z ~ 10) or the faint emission of 

conditions. Its 21-cm emission line is a key probe of the struc- clusters of galaxies at moderate redshifts (3 ^ z ^ 5) is shifted 

ture and dynamics of galaxies over a broad range of different ^ j^is soft X-ray band. It is not possible to analyze the X-ray 

^ densities and temperatures. H i is the key element to understand jata quantitatively without knowing in detail the distribution of 

- - the evolution of the ISM in general ( |Kalberla & Kerp| | 2009 ). 21- j^e Galactic interstellar medium across the target, 
cm line observations of the Milky Way are of particular interest 
because small scale structures are observable which can not be 

resolved with observations of any other galaxy. A number ofno- High resolution single dish data ar e needed but neither the 

table Items have been detailed by|McClure-Griffiths et al.| (120091 h i Parkes All Sky Survey (HIPASS) bv lBarnes etal](r2 00n nor 

hereafter Paper I), this includes the impact of massive stars on the reprocessed version (HIPASS-HVC) by'Putman et alj (i2002h 

the ISM, the disk-halo interaction, the ISM life-cycle and the provide accurate column densities for the Galactic emission in 

formation of cold clouds. H i is a perfect tracer of the Galactic the southern sky. The GASS provides such data. For the north- 

structure and there is a good correlation between gas and dust ^^ ^^y th e Galactic Arecibo L-Band Feed H i (GALFA-HI) sur- 

( |Schlegeletal. || 1998[) . vey (Peek_&Heiles|I^Qi) an d the Effels berg-Bonn Hi Survey 

Galactic H i data are important for an exploration of the early (EBHIS) are still in progress (Winkel et al. 2010). Major param- 

universe because of the photoelectric absorption caused by the eters for these surveys are compared in Table [II Currently the 

LAB and the GASS are the only surveys that are easily accessi- 

Send offprint requests to: P. M. W. Kalberla ble for data mining. 
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Parameter 



LAB 



HIPASS-HVC 



GASS 



GALFA-HI 



EBHIS 



full sky 6<2° 

36' 14.4' 

-450 < V < +400 -700 < v < +1000 

L27kms-i 26.4 km s"! 

90 mK 24 mK 



Sky coverage 
Average beam FWHM 
LSR velocity range 
Channel width 
IctTb noise 

Table 1. Survey parameters for major single dish surveys of the Galactic 21 -cm line emission. GALFA-HI and EBHIS are still in 
progress. 



-468 < V < +468 

l.Okms-i 

57 mK 



-r <(5<38° 

3.4' 

-750 < V < +750 

0.18kms-i 

< 100 mK 



(5>-5° 

9.5' 

-1000 <v< +1000 

1.25kms-i 

<90mK 



The Galactic All-Sky Survey (GASS) is a 21 -cm line sur- 
vey covering the southern sky for all declinations 6 ^ 1°. 
The observations were made with the multibeam system on 
the 64-m Parkes Radio Telescope. The intrinsic angular reso- 
lution of the data is 14^4 (FWHM). The velocity resolution is 
1.0 kms"^ and the useful bandpass covers a velocity range 
IvisrI ^ 468 kms"^ for all of the observations; some data cover 
up to IvisrI ^ 500 km s"^ . The GASS is the most sensitive, highest 
angular resolution large-scale survey of Galactic H i emission 
ever made in the southern sky. The first data release is available 
at |http://www. atnf . csiro . au/research/G AS S 

A detailed description of the GASS project, focusing on the 
survey goals and techniques has been given by Papier I, In Paper 
I all of the initial data reduction and imaging is described in 
detail. Here we focus on the post-processing, in particular on 
corrections for stray-radiation, instrumental baselines, and radio 
frequency interference (RFI). We will note where our data re- 
duction diff'ers from that described in .Paper I. 

For 21 -cm line observations there are a few unavoidable in- 
strumental issues that can potentially degrade the data. The first 
is that a radio telescope does not have a perfect beam, so some 
radiation enters through the sidelobes of the antenna, causing 
spurious emission features. The second eff'ect is that the receiver 
gain may drift, and it can suff'er from instrumental spectral struc- 
ture (the instrumental baseline) that may be time-dependent. 
These influences need to be minimized, but the raw observations 
do not allow a clear identification and separation of individual in- 
strumental problems. For example, spurious profile wings may 
be caused by stray-radiation but also by instrumental baseline 
defects. Temporal fluctuations in the line emission may be due 
to gain variations but may also be caused by spurious emission 
from the sidelobes. 

A strategy is needed to disentangle these instrumental is- 
sues. For the stray-radiation correction a basic prerequisite is 
the knowledge of the H i brightness temperature distribution 
provi ded by the L eiden/ Argentine/Bonn (LAB) 21 -cm line sur- 
vey (iKalberla et al. 2005); most important in the context of 
the GASS survey is the contribution from the 30-m dish of 
the Instit uto Argentino de R adioastronomia (lAR) in the south- 
ern sky (iBajaja et al.ll2005h . These surveys provide clean data 
with an angular resolution of ~ 36' and allow us to solve the 
Fredholm integral equation which describes the deconvolution 
of the raw observations (antenna temperatures) with respect to 
the influences of the antenna response. Bootstrapping from the 
LAB survey we use a model of the antenna response for the 13- 
beam receiver based on the known feed patterns and the tele- 
scope geometry. Section [2] describes how spurious features from 
the sidelobes are removed. 

Next we remove baseline defects. For this the LAB survey is 
used for initial estimates of velocity ranges that do not contain 
line emission. In Sect. [3] we describe how baselines are refined 
in an iterative way. A gain calibration requires some preliminary 



baseline corrections but such a calibration also aff'ects the ac- 
curacy of the corrections for stray-radiation. The calibration is 
described in Sect. IH All corrections discussed in Sects. [2] to |4] 
need to be iterated several times. 

RFI is increasingly a problem for sensitive observations, 
and post-processing mitigation strategies are essential. We use 
flags which were set in the first stage of the data processing by 
Livedata described in lPaper I and iteratively remove the remain- 
ing interference as discussed in Sect. [5] Throughout the reduc- 
tion we take care that the corrections described in Sects. [2l to |4] 
are not aff'ected by RFI. 

In comparison with the LAB survey, the GASS represents an 
improvement in velocity resolution and sensitivity, but most im- 
portant is the enhanced spatial resolution (Table [TJ. The genera- 
tion of FITS images with user-definable parameters is described 
in Sect. [6] We discuss some examples of the newly corrected 
data in Sect. [7] and complete the paper with a description of a 
web interface for download (Sect. (8]). 



2. The stray-radiation correction 

A major problem when observing 21 -cm emission lines is 
caused by the fact that Galactic emission is seen in all direc- 
tions. Most prominent is the bright Galactic plane which extends 
across the whole sky. It is unavoidable that stray 21 -cm emission 
from the plane will get picked up by sidelobes of the antenna. 
The main sidelobes are the result of reflections from the feed 
legs that carry the prime focus cabin. Also radiation originating 
from regions outside the rim of the reflector, th e spillover re- 
gion, is received. Stray-radiation varies with time (Ivan WoerdenI 
1962) and for positions with very low Hi the spurious emis- 
sion features_jnayainouiU_^^ 

nal (IKalberla et al.|[T98QaHLockman et al.lll986l:lHartmann et al.1 
Il996l) . In some GASS spectra stray-radiation was found to com- 
prise 35% of the observed emission. 

Stray-radiation may not be neglected if we are interested in 
low-intensity profile wings that may originate from lines with 
large velocity dispersion, or in deriving accurate Hi column 
densities from 21 -cm observations. Observed lines may origi- 
nate either from the main beam, the "real" signal, or may have 
been picked up from other directions, the "stray" signal. Without 
prior knowledge of instrumental issues it is also very difficult to 
distinguish between faint line emission and instrumental base- 
line defects. 

One way to minimize stray-radiation is by minimizing scat- 
tering surfaces in the telescope aperture. The optimal case ap- 
pears to_beaJelesco2econstruc^ with a clean unblocked aper- 
ture (iDickey &Lockmanlll990l). Prime examp les are the Bell 
Labs horn reflector antenna ("Stark et al." 1992") and the Green 
Bank Telescope (GBT; Prestage et al. 2009). A typical parabolic 
reflector antenna has a main beam efficiency of 70%, for un- 
blocked apertures this improves to ^ 90%. Sidelobe contribu- 
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tions are thus reduced from 3 0% to 10% but are still recognizable 
(iLockman & Co ndon! 1200 5*). Depending on the observed posi- 
tion, emission from the spillover region may cause severe con- 
taminations in the case of an u nblocked aperture w ith a receiver 
in the secondary focus (.Robishaw & Heilesll2009l) . 

A numer ical solution o f the stray-radiation problem was 
proposed by vanWoerden (1962). He demonstrated convinc- 
ingly that the stray component can be determined by convolv- 
ing the antenna pattern with the brightness temperature distribu- 
tion on the sky. Yet, neither the knowledge about sidelobe struc- 
tures, nor the computing power was sufficient enough for accu- 
rate calculations in 1962. The first reli able calculation s became 
available for the Eff'elsberg telescope ( Kalberla e t al.' 198 0a,b) 



and l ater for the resurfaced Dwingeloo dish (iHartman n et al.l 
1199 6*). the 26-m Telescope at the Dominion Radio Astrophysical 
Observatory (Higgs & Tapping 2000), and the 30-m telescope 
at Villa Elisa (Bajaja et al. 2005). A refined correction of the 
Leiden/Dwingeloo survey (Hartmann & Burton .1997 ) in com- 
bination with the Vill a Elisa observations led to the LAB survey 
(iKalberla et al.ll2005l) . Currently the LAB is considered to be the 
most accurate Galactic all-sky survey. 

2.1. Basics 



We use a correction algorithm developed by Kalberla et al.l 
(1980a). Antenna temperatures Ta observed by radio telescopes 
can be approximated in Cartesian coordinates as a convolution 
of the true temperature distribution T on the sky with the beam 
pattern P of the antenna 



r,(,,„ = /p<,-/,,-/,r,/,/M/./^ 



(1) 



This is an approximation, but is sufficient to demonstrate the ba- 
sics. In general, Ta is time and frequency dependent, spherical 
coordinates should be used and the integration needs to be ex- 
tended over the observable part of the sky as well as the ground 
if there are reflections. We do not include in these equations the 
effects of atmospheric attenuation, though this must be consid- 
ered for emission entering both the main beam and the sidelobes 
(ILockm an 2002). 

For the pattern P of the antenna we use the normalization 



/ 



P(x, y)dxdy = 1 



(2) 



in addition, we split the pattern into the main beam area (MB) 
and the stray-pattern (SP) 



r,(x,j)= r p{x-x\y-y)T{x\y)dx'dy + 

J MB 



JSP 



P(x-x\y-y)T(x\y)dx'dy' 



(3) 



rjMB = I P(x,y)dxdy 
Jmb 



Defining the main beam efficiency tjmb of the telescope as 

(4) 

we may rewrite Eq. 3 as 
TB(x,y)=^^^^^-— [ P(x-x\y-y)T(x\y)dx'dy.(5) 

r/MB riMB JSP 

Replacing Tb by T leads to a Fredholm integral equa- 
tion of the second kind which can be solved for riMB > 0.5 



(iKalberla et al .1 1 1 9 SOah . Tb is the so-called brightness tempera- 
ture, resulting from a convolution between the main beam of the 
telescope and the true sky temperature T. A solution of Eq.[5]is 
simplified considerably for a known all-sky brightness tempera- 
ture distribution T. We use profiles from the LAB survey as the 
currently best possible approximation to T. For a solution of Eq. 
[51 P needs to be determined. This is described in Sects. l23l and 
[241 



2.2. The correction algorithm 

The nii merical correction algorithm is based on IKalberla et al.l 
(Il980ah . It was originally developed for the Eff'elsberg 100-m 
telescope, wi th some improveme nts it was then applied to the 
LAB survey ("Kalb erla et al.ll2005l) . We modified this procedure 
for use with multi-beam systems, allowing individual correc- 
tions for each of the beams. Due to the rotation of the off'set 
receivers with respect to the telescope structure, we take into ac- 
count the variation of each receiver's beam pattern on the sky 
with time. To improve the accuracy of the calculations we also 
changed the convolution algorithm for the inner part of the an- 
tenna response pattern (within 6° of the main beam). From the 
LAB survey, a new and more accurate input sky for the convolu- 
tion (Eq.O was generated. 

2.3. The multi-beam antenna pattern - near sidelobes 

To solve Eq. [5] it is necessary to know the complete antenna 
response P for each of the 13 beams of the Parkes multi-feed 
system. In practice there are severe limitations in measuring the 
sidelobes to the required accuracy. We have therefore chosen to 
model the sidelobes. In general, the far-field pattern of an an- 
tenna is the autocorrelation of the aperture plane distribution. 
Using the measured feed response pattern and the telescope ge- 
ometry, including shadowing caused by the focus cabin and the 
feed support legs, we derive the complex aperture distribution 
function for each feed. The pattern can then be calculated by 
Fourier transformation. 

For the 13 beams we distinguished three characteristic cases, 
the central beam and one feed from the inner and outer ring of 
feeds (Fig.[T]top). The feed horns are equipped with two pairs of 
orthogonal-linear probes, which, while in parallactic tracking, 
rotate relative to the feed support legs. In turn the complex aper- 
ture distribution function changes with time. We approximate 
this rotation with an accuracy of ±5° . The fact that the three feed 
support legs cause a global six-fold symmetry for the antenna 
pattern simplifies the situation considerably. We need only to 
calculate 5 individual patterns for each of the rings with radial 
feed off'sets of 29f 1 and 50f8, respectively. 

Fig. [T] shows the feed configuration for the multibeam sys- 
ten£l and displays the antenna response pattern within a radius 
of 6° for feeds 1, 2, and 8. We see characteristic changes due 
to the beam displacements. The main beam develops a coma 
lobe that increases in strength with feed off'set. Sidelobe struc- 
tures become asymmetric and have a more pronounced ringing 
with increasing feed off'set. An important feature of the multi- 
feed system is that the azimuthal asymmetries in the antenna 
pattern due to the feed off'sets decrease on average with increas- 
ing distance from the main beam. Extrapolating this property 
to far sidelobes, off'set from the main beam by tens of degrees, 
implies that diff'erences between individual feeds should vanish. 
Indeed we verified this behavior as described in Sect. I2.4[ 



http://www.atnf.csiro.aU/research/multibeam/.overview.html 
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Fig. 1. Feed configuration for the receiver (top) and antenna pat- 
terns within a radius of 6° for beam 1, 2, and 8 (from top to bot- 
tom). The display is restricted to levels between -20dB (white) 
to -50db (black). 



In the correction algorithm it is not necessary to use a pattern 
with all the details shown in Fig. [T] as the LAB data were mea- 



sured on a grid of only 0?5 x 0?5. We have therefore modeled the 
inner 6?0 of the antenna patterns on a grid containing 468 cells 
in cylindrical coordinates. The cells have an azimuthal extent of 
10° and a radial extent of 0?2 out to 2? 2 from the main beam, a 
radial extent of 0?3 out to 2? 5, and of 0?5 thereafter. 



2.4. The far sidelobes 

The far sidelobes arise predominantly by reflections from the 
feed support legs, causing the so-called stray-cones, and at the 
rim of the telescope, the spillover. The location of these struc- 
tures on the sky was determined from the telescope geometry 
and verified by numerous tests made with a transmitter. The side- 
lobe levels from the spillover lobes were estimated from the edge 
taper of the receiver feed. For the stray-cones we used estimates 
from the the aperture blocking by the feed legs. These estimates 
were also based on sidelobe levels found previousl y for the tele 



scopes at Eff^elsberg, Dwinge l oo and Villa Elisa (Kalb erla et al 



1980at lHartmannetaDll996l: iBaiaiaet al.l l20Q5: Kalb erla et al 



2005). The far sidelobe components were adjusted individu- 
ally in a way similar to that described by Kalberla et al. (2005h . 
searching for a consistent solution of Eq. [S] 

2.5. Ground reflections 

Sidelobes that touch the ground receive broad-band thermal 
noise, but to first order no line radiation is expected. However, 
depending on the ground cover, reflections may happen. While 
scrub acts as an absorber, low grass and soil can reflect 21 -cm 
line radiation. We determine this contribution by surveying the 
ground around the telescope, and then using ray-tracing methods 
to estimate the efl'ective sidelobe response. We used an albedo of 
0.2. 

3. Removal of instrumental baselines 

3.1. Preliminary considerations 

In many cases, stray-radiation produces extended profile wings 
that mistakenly may be interpreted as instrumental baseline 
structure. A correction for stray-radiation could cause negative 
intensities if applied after the baseline correction. It is there- 
fore essential to apply baseline corrections only after the stray- 
radiation is eliminated. For this reason we had to deviate from 
the processing described in iPaper ll : for the first data release no 
correction for stray-radiation was applied. We used Livedata as 
described in their Sect. 2.3 for the bandpass correction, but with- 
out the post processing for baselines according to Sect. 2.3.3 in 
iPaper it The algorithm described below is in lieu of the post- 
processing described in lPaperJ . 

After scaling the reduced data to antenna temperatures cor- 
rected for atmospheric absorption (Sect. |4]) the stray-radiation 
was subtracted and afterwards the instrumental baseline was cor- 
rected. Flags set by Livedata or our own software to indicate the 
presence of RFI were taken into account, and aff'ected channels 
were disregarded for all of the baseline fits described below. 

3.2. The algorithm 

For each spectrum of the stray-radiation-corrected GASS 
database we used the LAB survey for an initial estimate of the 
velocity extent of the expected Hi emission. After smooth- 
ing the interpolated LAB profile to an eff'ective resolution of 
Avi/2 = 5 kms"\ we used velocities where T^ab < 0.9 K as 
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the range for an initial fit of the instrumental GASS baseline. In 
addition we folded the LAB profile in frequency to mimic the 
second sp ectrum used by Livedata to determine the instrumental 
bandpass (iPaper \ Fig. 3). Flags were transferred accordingly. 

GASS spectra from both polarizations were averaged and 
smoothed to Avi/2 = 5 kms"^ for each of the individual 5 sec 
dumps. Using the initial estimate for the emission-free part of the 
profile we determined the instrumental baseline by fitting a 9^^ 
order polynomial. Any channel with signal exceeding four times 
the rms noise over the fitted region was considered to contain 
possible H i emission and was then excluded from the baseline 
region. The baseline fit was repeated with an 11^^ order poly- 
nomial, again searching for potential emission that needed to be 
excluded from the fit. As a measure of the fit quality we deter- 
mined the residual rms scatter within the baseline region. For 
a few percent of the data we found that a 9^^ order polynomial 
provided a better fit. Accordingly we used this fit. Finally we 
took into account that weak profile wings with intensities below 
the automatic Act cut-off' level could aff'ect the fit parameters. To 
avoid such a bias we excluded 6 more channels on each side of 
the H I emission from the baseline region that was fitted. 

After these initial steps we considered both spectra for each 
polarization individually. The spectra were smoothed and fitted 
with the parameters determined previously to find out whether 
some more features needed to be excluded from the baseline re- 
gion. After that the final baseline fit was applied to the original 
observations. 

The algorithm and the parameters described above were 
tested and optimized by varying the polynomial order and try- 
ing other functional forms. We particularly tried to avoid high 
order polynomials. As described in Paper I, the average quo- 
tient spectrum for each of the scans (visible as patches in Fig. 
P botto m) was fitted with a 15^^ order polynomial. Figure 3 of 
iPaper J demonstrates that emission features are not aff'ected by 
this fit. After this global bandpass fit there remain some time de- 
pendent ffuctuations that need to be corrected independently for 
each spectrum with 5 sec integration time. A first guess would be 
that a similar polynomial is needed to get rid of such ffuctuations. 
This is fortunately not the case. In the initial Livedata reduction 
for the ffrst data release a 10^^ order fft for post-processing of in- 
dividual spectra was sufficient. Our results conffrm this ffnding 
with minor modiff cations for our procedure. We tested lower or- 
der polynomials. These led for many regions to obvious baseline 
defects. We found spurious features in emission free regions that 
could only be avoided by allowing higher orders. 

Assuming that some parts of the residual baseline ripples 
could be caused by standing waves, we also tested sine wave 
ffts in addition to polynomials. Sinusoid solutions were found 
to be insigniffcant and were discarded. We carefully checked 
whether errors in the LAB data might aff'ect the GASS baseline 
correction, but the iterative reffnements that were applied suc- 
cessfully avoided any biases. We tested also whether the high 
order polynomial ffts might have aff'ected weak emission fea- 
tures. We found no indications for such problems. In particu- 
lar the RFI feature observed in March 2006, appearing similar 
to a HVC (Sect. 15.5b , remained unaff'ected by the fft. Also test- 
ing whether proffle wings in data that remained uncorrected for 
stray-radiation could be removed by baseline fftting did not in- 
dicate an unwanted overcorrection by the baseline fft (see Sect. 

m. 

Some residual systematical baseline problems remain in a 
few regions, on average typically at a level of ^ 40 mK but oc- 
casionally up to 100 mK. This is close to the rms uncertainties 
of ~ 60 mK for individual maps as discussed in Sect. [71 We were 



Table 2. Tb calibration positions and parameters 



Pos. 



1 



Peak Tb 

(K) 



Area 
(Kkms- 



Vel. range 
(kms-^ 



S8 


207? 00 


-15?00 


75.6 


879.5 


-5.1 to +22.3 


S9 


356?00 


-4!00 


83.5 


930.8 


-1.5 to +15.2 


S6 


1?91 


41?42 


51.0 


285.5 


-5.8 to +4.8 



unable to identify the reason for such systematic baseline errors. 
They may be caused by RFI or bad ffts but also by uncertainties 
in the stray-radiation correction. 



4. Calibration 

An accurate correction for stray-radiation demands a careful cal- 
ibration. Errors in the antenna temperature scale cause a mis- 
match with the stray-radiation correction term in Eq.[5j The side- 
lobe response is frequency dependent; calibration errors, there- 
fore, cause not only an increased scale error but also aff'ect the 
proffle shape. We used regular observ ations at the standard po- 
sitions SB, S9, and S6 (IWilliamslll97 3) and determined average 
calibration factors for each of the observing sessions. An ini- 
tial temporary calibration was determined from uncorrected data 
(Paper I). After correcting the data for stray-radiation we rede- 
termined the Tb calibration factors and applied these new factors 
for the following calculations. 

Table [2] lists the ffnal calibration data; the brightness tem- 
perature scale is matched to a common temperature scale ffrst 
established by lWilliaini (Il973l) . The brightness temperature for 
most of the major surveys are matc hed through this reference ; 
this includes the Hat C reek surveys (Weaver & Williams "1973'; 
iHe iles & Habing|[T974l) . the Le iden-Green Bank survey (Bu rtonI 
19851), southern sky surve ys bv lKerretall (Il98lh . the BeU Labs 
survey ("Stark etal." 1992"), the LAB survey, and Eff'elsberg and 
Arecibo data (Peek & Heiles 2008; Winkel et al. 2010). 

A few calibrati on regions have been mapped by 
iKalberla et al.l (Il982h . They calculated the beam response 
in detail and provided also calibration factors as a function of 
the beam shape. We use factors appropriate for a 14.4'beam. 

The easiest way to determine the calibration factors is to 
use peak temperatures, bi it line integrals over a limited velocity 
range are more accurate (Williams" 1973) and we preferred this 
method. Please note that velocities listed in Table [2l are integra- 
tion limits while velocities given by Williams (19 731) ar e center 
velocities for the edge channels (see Kalb erla et al.i 1982 1). Figure 
[21 (left) compares overall gain factors applied to the initial cali- 
brations (blue) with those after completing the ffnal corrections 
(red and green). The overall calibration uncertainties (right) for 
most of the feeds could be reduced to typically 1.7%. Only two 
receivers were less stable with an rms scatter of 3% and 3.5%. 
The uncertainties are of statistical nature. Combining observa- 
tions with several beams at diff'erent seasons we expect average 
temperature scale uncertainties well below 1 % in the ffnal data. 
We note that interference can cause a loss of calibration in some 
spectra that would be difficult to detect as it would appear as a 
scale factor over an entire spectrum. We believe that this eff'ect 
is small in the overall survey but may be present in some ffnal 
spectra. 
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Fig. 2. Left: Gain factors applied to the initial calibration (blue stars and purple squares) and cross check after the final data reduction 
(red plus signs and green X's). Right: RMS scatter in gain calibration for individual feeds, before and after correction for stray- 
radiation (same symbols as in left panel). 



5. RFI mitigation 

5.1. Replacement of data flagged by Livedata 

Narrow-line interferenc e was flagged during the basic data re- 
duction with Livedata (Paper I, Sect. 2.3.2). Most of this type 
of RFI occurs at a fixed topocentric frequency, so as the tele- 
scope scans across the sky the RFI moves in successive channels 
in the LSRK frame. We use the kinematical (radio) standard of 
rest definition with a velocity of 20 km s"^ in direction ra,dec = 
[270°, +30°] (epoch 1900.0). Most of the features are observed 
to stretch over large angular distances and aff'ect fields that have 
been observed independently at diff'erent seasons (Fig. O. This 
behavior implies that the RFI remained approximately constant 
for long periods in time, aff'ecting most of the data. 

Profiles with more then 30 flagged channels were completely 
discarded during imaging; for less serious contaminations the 
flagged data were replaced. After correcting profiles for stray- 
radiation and instrumental baseline, we interpolated LAB sur- 
vey data to replace the flagged channels. Alternatively we have 
tried to replace data by interpolating the GASS spectra using 
the nearest ten good data points. Using LAB data resulted in all 
cases in much cleaner maps. Still, some weak stripy features re- 
main that appear to move across the sky in successive velocity 
channels but mostly these artifacts remain within the noise (Fig. 
[3]). About 0.1% of all data have been flagged by Livedata and 
were replaced by interpolated data in this first step. We kept the 
Livedata flags to avoid such interpolated data getting propagated 
by the RFI mitigation described in the following sections. 

5.2. Median filter RFI rejection 

Even after removal of the channels flagged in the initial data re- 
duction, numerous instances of RFI remain in the data. Some ap- 
pear as point sources, but the most obvious show "footprints" of 
the 13-beam system that move on the sky in scanning direction 
(Fig|4]). As before, the motion in successive velocity channels 
is caused by changing LSR projections of the fixed Topocentric 
RFI frequencies. Features of this kind were found previously in 
the HIPASS and are described in detail bv 'Barnes et al.l (l200ll 
Fig. 7). For the first data release the remaining RFI was removed 
by median gridding when calculating FITS images (JPaper 1 
Sect. 2.3.5). 



Here we follow a diff'erent strategy, taking advantage of the 
fact that in GASS typically 40 individual spectra contribute to 
the final data in every resolution element. First, for any observed 
position in the sky we selected all observed profiles within an 
angle e of typically a few arc-minutes. From these data we cal- 
culate the mean for each channel, the corresponding standard 
deviation, and the median, disregarding data flagged either by 
Livedata or our own post-processing. The average rms noise 
level Trms = o-{T) and its one sigma deviation (t{Trms ) from 
the average was determined for the low level emission part of 
the profile {Tb < 0.5K). 

We have two tests for RFI. In the first, channels with rms 
deviations 1ST rms > ^o-(Trms ) were considered as highly suspi- 
cious for RFI contamination and investigated in more detail. We 
use as a second indication of RFI circumstances where the mean 
and median diff'er from each other by more than one standard 
deviation. 

To eliminate outliers for the channels that were selected in 
both of the two tests we determined if individually observed 
brightness temperatures deviate from the median by more than 
two times the standard deviation. In the case of a purely ran- 
dom distribution this would exclude 5% of the good data, and 
we consider this an acceptable cost to eliminate true RFI. We 
flag outliers and repeat, calculating new mean and median esti- 
mators and once more excluding any outliers. 

Our two sigma limit may be compared with the more rig- 
orous method for the elimination of suspect data proposed by 
iPeircd (Il852h which is based on probability theory. Accordingly 
a rejection of a data point that deviates from the mean at a two 
sigma level is justified if there are 13 independent observed data 
values. For 6 = 6' (see below) we have on average 40 individ- 
ual profiles and a two sigma limit is appropriate if about 10% 
of the data are suspect. This condition is valid for many cases 
but we find also situations where 20% of the measurements are 
outliers. In this case a cutoff at 1.6 (t (Trms) would be more ap- 
propriate according to iPeircd (Il852h . We deviate from Peirce's 
criterion in two ways; to simplify the calculations we use a fixed 
2 o-(Trms ) cutoff' and we define outliers by their deviations from 
the median and not, as originally proposed by Peirce, from the 
mean. The latter is essential since outliers caused by RFI can be 
found at ten to hundred times the standard deviation, a circum- 
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Fig. 3. Typical example for extended interference patterns that were flagged by L ivedata in the first stage of the data reduction. Left: 
as observed, right: after RFI removal using LAB data as described in Sect. 15 11 The intensity scale is linear for -0.2 < T < 0.5 K. 
Some weak structures remain but are mostly at |r| < 40 mK 



stance that was not considered by Peirce; in such cases the mean 
is usually seriously aff'ected. 

To replace values outside the two sigma range, the best avail- 
able estimate is the median value scaled by the proper beam 
weighti ng function, for deta ils we refer to the extended discus- 
sion by lBarnes et al.l (l200ll Sect. 4.2.2). We replace the previ- 
ously selected data points by the weighted median and set a 
flag to distinguish medians from observed values. A fraction 
~ 5 • 10"^ of the data is afl'ected by RFI and replaced by the 
weighted median. 

One of the criteria that we used to detect RFI was the fact 
that mean and median should not difl'er by more than one stan- 
dard deviation. We applied this criterion also to test whether the 
median estimator is consistent with the mean of all data after 
exclusion of the outliers. Both methods were found to produce 
nearly identical results. An example of RFI in scanning direction 
and the result after application of our median filer is given in Fig. 

n 

Figure [5] shows spectral details for one position in this field. 
Displayed are the mean (red), the rms scatter (green), and the 
median (blue) derived for all profiles that are closer than e < 6' 
to the position / = 0°, /? = -13?5. Data in the left plot suff'er from 
strong interference at visr ~ 157 kms"^ 55 ^ visr ^ 58 kms"^ 
and at visr ^ 400 kms"^ Only a part of the data was flagged 
by Livedata and replaced by values from the LAB. Despite all 
defects the median (blue) is essentially unafl'ected by RFI. The 
data after RFI elimination are shown in the plot at the right side. 
Profiles for mean and median are within the noise identical, the 
rms is close to the expected value of 0.4 K and shows no strong 
enhancements. We find some fluctuations close to the line emis- 
sion at Visr ~ kms"^ This may have been caused by genuine 
fluctuations in the H i line emission. To avoid a possible degra- 
dation of such fluctuations no automatic RFI cleaning was ap- 
plied for Tb ^ 0.5K. These two plots demonstrate how greatly 
the rms is afl'ected by RFI, and how deviations between median 
and mean are obvious. Our RFI filter is triggered by such devia- 
tions. 



5.3. The beam-specific parameter e 

A critical parameter for RFI filtering is the radius e for a com- 
parison of the data at adjacent positions. A large value provides 
a large statistical sample but at the same time point sources or 
steep gradients in the brightness distribution may be aff'ected by 
the filter. We found in tests that 6 = 6' is sufficiently robust. 
Tests with H i point- sources have shown that larger values could 



degrade t he source distributio n. Our result is in excellent agree- 
ment with lBarnes et al.l (1200 ih . and our median estimator is iden- 
tical with the parameter Vmax = ^' that was used for HIPASS. 



5.4. RFI in emission line regions 

RFI that falls in channels containing H i emission was usually 
removed successfully by the above procedure, except in regions 
with strong spectral gradients. Examples of such critical cases 
are positions with strong absorption lines. Sources with a contin- 
uum flux exceeding 200 mJy were excluded from RFI rejection. 

To avoid any possible degradation of the signal due to an au- 
tomatic RFI filtering, the median filters described in Sect. 15.21 
were restricted to data with low brightness temperatures. We 
used a limit of T^ax ~ 0.5K. The remaining data were checked 
visually. Only a few cases showed some residual RFI. For these 
we increased the filter limits such that only channels with suffi- 
cient large rms deviation ATrms were aff'ected. We verified that 
the remaining data were unaff'ected by the median filter. As a 
general rule we tried to avoid as far as possible any modification 
of the H I emission lines through the RFI filtering process. This 
implies, however, that some RFI may remain undetected and be 
present in the final data. 

Some of the spectra were aff'ected so severely by RFI that 
they needed to be removed completely. We eliminated all pro- 
files whose mean rms within the baseline region exceeded the 
average noise level at its position by a factor of three. About 
0.3% of the observed profiles were rejected for this reason. 



5.5. RFI in March 2006 

For a short period between 18 and 21 March 2006 observations 
were degraded by broadband RFI with quite diff'erent charac- 
teristics than any of the narrow band RFI described above. The 
weak spurious lines had approximately Gaussian shapes with 
center velocities in the range 270 to 3 10 km s"^ and velocity dis- 
persions of about 15 km s"^ . None of the RFI rejection methods 
described above applied to this case. These spurious lines were 
identified by fitting Gaussians; all aff'ected channels were then 
flagged. Baseline corrections were repeated, finally all flagged 
data were replaced by medians as described before. This way 
we recovered data that were discarded in the first data release. 
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Fig. 4. HI emission at visr = 61.8 kms"^. Data heavily affected by RFI in scanning direction of the telescope (left) are compared 
with the clean map (right). The intensity scale for the range -.04 < T < 20 K is logarithmic, emphasizing low level emission. 
Yellow colors correspond to a brightness temperature of ~ IK. 
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Fig. 5. Data at the position / = 0°, /? = -13°5 before and after RFI mitigation. The average is plotted in red, the rms scatter in green, 
and the median in blue. 



5.6. Removal of bandpass ghosts 

To maximize sensitivity and recover all the extended H i emis- 
sion, the GASS was made using in-band frequency- switching 
with an offset of 3.125 MHz corresponding to 660 kms"^ The 
Livedata bandpass correction causes for every real emission 
line feature an associated, spurious, negative image ("ghost") in 
the other band, displaced by ±660 kms"^ Most of these nega- 
tive images fall outside the velocity coverage of the survey but 
high velocity clouds (HVCs) with emission lines at |visr| ^160 
km s"^ cause ghosts. In particular, the Magellanic System causes 
a strong and extended bandpass ghost at negative velocities. 

As Fig. 3 of iPaper J shows, ghosts appear only in one of 
the two IFs; the other band is unaffected. It is possible to pro- 
duce maps without ghost features if one uses only one IF for 
IvisrI ^ 160 kms"^ A second possibility is to flag ghosts and 
eliminate them with a treatment analogous to one of the the RFI 
mitigation strategies discussed above. Accordingly the ghosts 
were replaced by medians determined from the alternate ghost- 
free IF band. Such maps have the advantage that they provide a 



better sensitivity for most of the field but remaining biases can 
not be excluded. We determined for both bands all HVC emis- 
sion features with |visr| ^ 160 kms"^ and Tb > 0.2 K (corre- 
sponding to 3cr(TRMS ) for 6 = 6') and flagged ghosts in the other 
band. These features were treated analogously to the broadband 
RFI discussed in the previous section. 

6. Imaging 

For HIR\SS, as well as for the first data release of the GASS, 
imaging was performed by Gridzilla using weighted median es- 
timation which is robust against RFI and other artifacts. For the 
second data release, we use separate and independent procedures 
for RFI rejection and imaging. As described in the previous sec- 
tion, our median filter is largely consistent with the HIR\SS me- 
dian filter algorithm for a radius 6 = 6'. However, the major dif- 
ference is that we apply this filter only to individual channels if 
two conditions apply: 1) an exceedingly large scatter is observed 
with a 3 times larger dispersion than the mean dispersion caused 
by thermal noise and 2) mean and median differ by more than 
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Fig. 6. Effective FWHM resolution W of GASS images as a 
function of the user-definable Gaussian interpolation kernel Wk. 
For of Wk > 7' we find W ~ W9 + 0.666W)^. For comparison we 
show the eff'ective resolution obtained by median gridding and 
by Bessel interpolation (horizontal lines). 
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Fig. 7. Eff'ective rms uncertainties of GASS channel maps with 
1 km s"^ velocity resolution as a function of the user-definable 
Gaussian interpolation kernel. Crosses indicate fitted values de- 
termined at high Galactic latitudes, the solid black line repre- 
sents upper limits of the noise, valid for most of the positions in 
the Galactic plane. For comparison we show the eff'ective noise 
at high latitudes obtained by median gridding and by Bessel in- 
terpolation (horizontal lines). 



one standard deviation. Only data that deviate from the median 
by /ST RMS > 2 • o-{Trms) are replaced by the median estimate. 

For the data release described here the gridding was per- 
formed using a two-dimensional Gaussian with a user-definable 
dispersion. For computational reasons we disregard data out- 
side a three sigma cutoff' radius. The interpolated value for that 
pixel is the weighted mean over all available input data. Except 
for the irregular distribution on the sky of the input data, this 
interpolation is equivalent to a linear convolution of the ob- 
served brightness temperature distribution with the kernel func- 
tion of the gridder. The resulting eff'ective beam is approximately 
Gaussian. For a source width Ws, an eff'ective average telescope 
beam width of Wt, and a width Wk of the user defined kernel we 
expect a resulting FWHM width 



We tested this relation by generating a series of images 
of the ultra-compact high- velocity clo ud H VC289+33+251 
(Putman et al. 2002; Briins & Westmeier 20041) . We used a grid 
spacing of 4f8, corresponding to 1/3 of the FWHM beam size. 
This object is isolated and has a FWHM diameter of Ws = 4 '4. 
Since it is unresolved by the Parkes telescope at 1.4 GHz the 
eff'ective beam W can easily be fitted. This was done for 18 
FITS cubes with 3' < Wk < 20' at two velocities, v = 250.7 
kms"^ and v = 251.5 km s~ ^ We deter mined first an eff'ective 

telescope beam width Wt = Jw^ - W^ = 14:2 ± i:0 from the 

known width of kernel and source. This value agrees well with 
the expected average FWHM beam width. The individual beams 
have 14' < FWHM < 14:5, with a mean of Wav = 14:4. 

Next we deconvolve the observed and fitted FWHM width 
Wobs of the image for the finite width Ws = 4:4 of the source, 

W = Jw^^^ ~ ^?- III Fig- we plot W as a function of the 
convolving kernel width Wk. To 3. good approximation we find 
for Wk ^ 6' the simple relation W = 10:9 + 0.666 • Wk. For a 
kernel Wk = 8' we obtain W = 16', identical with the resolution 
of the imag es calcii lated from weighted median estimators as 
presented in lPaper J . 

The noise level of the output map depends on the user de- 
fined smoothing kernel. We measured the noise in the same 
FITS cubes that we used to fit beam widths, at a velocity of 240 
km s"^ where neither H i emission nor RFI is visible. Figure [7] 
gives the results. The noise strongly depends on the smoothing 
kernel; for Wk = 7:0 we obtain an rms brightness temperature 
noise of 57 mK, comparable to the results in Paper I. The corre- 
sponding eff'ective resolution is W = 15:6. 

The noise level depends on the strength of the line emission. 
Within the Galactic plane for |Z7| ^ 1° we find a typical increase 
by a factor of 1.75. This approximate upper limit is plotted in 
Fig. [71 For a detailed map of the uncertainties we refer to the 
bottom map in Fig.|9l In this case Wk = 12:0 was chosen, result- 
ing in a typical noise level of 35 mK. It is obvious that the noise 
in the final maps is aff'ected by RFI mitigation and observational 
setup. Overlap regions for the raster chosen are clearly visible. 

We tested also the optimal gridding of single dis h on-the- 
ffy observations as proposed by IMangum et al.l (l2007l) with the 
function 



Ji(r/a) 2 

— — — exp-(r//?) . 
(r/a) 



(7) 



W 



^[m 



+ Wf + wl 



(6) 



/i is the Bessel function, a = 1.55, and b = 2.52 for a grid 
of 1/3 FWHM . This is the default for the AIPS task SDGRD 
(iGreisenll 19981) . The resulting resolution and noise are plotted as 
horiz ontal line s in Figs. [6l and [T] Both median gridding accord- 
ing to 'Paper T and optimal Bessel filtering are consistent with a 
Gaussian gridding for a kernel of Wy^ ~ 7:5. 



7. Quality of the Data 

To aid in the discussion of this release of the GASS data we 
have chosen a channel map at visr = -30.5 kms"^ (Fig. [8]). 
At this velocity stray-radiation eff'ects are strong and residual 
problems easily visible. In the top panel of Fig. [H we show a 
map that was generated without a stray-radiation correction but 
with the instrumental baselines removed as described in Sect. 
O In the middle panel, we plot the map from the final GASS 
database, corrected for stray-radiation and instrumental baseline. 
We have chosen a FWHM kernel Wk = 12:, resulting in rms 
uncertainties of ~ 35 mK at intermediate Galactic latitudes. In 
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Fig. 8. HI emission at visr = -30.5 kms"^. The intensity scale for the range -0.12 < T < 50 K is logarithmic, emphasizing low 
level emission. Top: GASS, corrected for instrumental baseline but without stray-radiation correction, middle: GASS corrected for 
stray-radiation and instrumental baseline, bottom: LAB survey for comparison, showing faint patchy spurious features that do not 
appear in the new GASS data, e.g. at / = 335°, b = -40°. 
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Fig. 9. Top: Stray-radiation contribution removed from the GASS at visr = -30.5 kms ^ the intensity scale is identical to Fig. [HI 
Bottom: rms noise in the final GASS data, corresponding to the intensity map shown in the central panel of Fig. [H 



the bottom panel of Fig. [HI we show the same channel from the 
LAB database for comparison. Figure [9] (top) contains the stray- 
radiation contribution removed from the middle map in Fig. [H] 

The bottom panel of Fig. [9] displays a compilation of the 
residual Icr rms uncertainties for a FITS cube generated with 
FWHM kernel Wk = 12f These uncertainties depend on the inte- 
gration time and on the system noise. Enhanced noise is caused 
by low elevation observations, by continuum emission from the 
Galactic plane, or by RFI mitigation (rejecting data). 

The final GASS map. Fig. [H] middle, shows no obvious ev- 
idence of residual stray-radiation problems. Time dependence 
of stray-radiation would cause a correlation with the observa- 
tional grid that is easily visible in Fig. [9] (bottom). Such resid- 
ual problems can be seen in the LAB channel map (Fig. [H] bot- 
tom). Individual blocks, corresponding to a grid of 5 x 5 po- 
sitions, show up at latitudes of b ~ -50°, a problem already 
mentioned by JKalberla et al.l (l2005l) . Similar block structures ex- 
isted also in the first version of the Leiden/Dwingeloo Survey 
(iHartmann & Burton 1997; H artmann et al. 1996 ) but could be 
removed by an improved stra v-radiation correction in the sec- 
ond data release of the LAB (JKalberla et al.ll2005b . Time vari- 



ability is caused predominantly by a seasonal variation of the 
stray-radiation and n eeds to be minimized by the correction 
(iKalberla et al.lll980allbh . The final GASS database is much less 
aff'ected by residual instrumental problems than the LAB survey. 
This implies in particular that there was no error propagation 
from the LAB to the GASS. We conclude that our corrections 
successfully removed most of the sidelobe influences, remaining 
systematical errors are probably below 1%. The performance of 
the survey is equivalent to observations with a telescope having 
^99% main beam efficiency. Scale uncertainties, caused by gain 
fluctuations, are also probably below 1% (Sect. [4]). 

The top panel of Fig. [H] shows that stray-radiation can not be 
removed with only our baseline correction algorithm as most of 
the stray-radiation remains in the map. This also demonstrates 
that our baselining algorithm is safe in the sense that it does not 
noticeably afl'ect extended profile wings. 

We searched for remaining systematical uncertainties in the 
instrumental baseline. Some residual uncertainties caused by 
RFI are visible in individual channel maps. In most cases such 
defects remain on average below 40 mK but there are a few more 
serious cases. In critical cases we recommend comparing GASS 
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data with the LAB database. Also a comparison of the cleaned 
database with the dirty data, offered by the web interface, may 
be helpful. 

8. Data products 

Data from the first release are accessible at 
|http ://w w w. atnf . csiro . au/research/G AS S/Data.html We provide 
a web interface at http://www.astro.uni-bonn.de/hisurvey for 
the second release to retrieve spectra and column densities 
for individual positions or complete FITS cubes. The user 
can specify Galactic or Equatorial coordinates, the velocity 
range, and the FWHM of the smoothing kernel for the optimal 
resolution. Default values for the gridding are provided as well. 
Some restrictions to the sizes of the FITS cubes apply, and we 
currently provide maps only in Cartesian projection. 

The user can check the data for residual contaminations by 
RFI or stray-radiation by downloading FITS cubes with dirty 
data (corrected for baseline and stray-radiation but no RFI filter- 
ing applied) or stray-radiation corrections that were subtracted 
as described in Sect. [2j Data at high velocities |visr| ^ 160 
kms"^ can be generated alternatively from both IFs or from 
single "ghost- free" IFs only (see Sect. 15.61) . In the latter case the 
rms noise for channels at velocities |visr| ^160 km s"^ increases 
by a factor of ^|2. 

Despite all eff'orts some instrumental problems may still exist 
in the second data release. Our web based policy for the on-the- 
fly generation of FITS cubes allows appropriate corrections. We 
ask users of the GASS data to report problems. Information on 
updates of the database will be provided on the web page. 



9. Summary 

The Parkes Galactic All- Sky Survey (GASS) measured the 
Hi emission in the Southern sky with complete sampling at 
all declinations ^ < 1° using the Parkes Radioe Telescope. The 
survey has an eff'ective angular resolution of 14^4 at a velocity 
resolution of 1.0 kms"\ and the data were obtained in a way 
that is sensitive to extended diff'use emission as well as compact 
sources. The GASS database contains 2.8-10^ individual spectra 
with 5 second integration time observed by on-the-fly mapping. 
The first data release together with a det ailed de scription of the 
survey goals and techniques was given in lPaper J . Here we focus 
on the post processing to correct instrumental eff'ects including 
stray-radiation, radio frequency interference, and residual instru- 
mental spectral baselines. 

The GASS is the most sensitive, highest angular resolution 
survey of Galactic H i emission ever made of the Southern sky. 
After correcting for instrumental eff'ects and RFI the GASS data 
are in excellent agreement with the LAB survey. This is impres- 
sively demonstrated with Fig.fSj FITS data cubes are available on 
request from http://www.astro.uni-bonn.de/hisurvey, Our soft- 
ware allows a flexible generation of FITS datacubes with efl'ec- 
tive resolutions 16 ^ FWHM ^ 25' (Fig.[6l). The correspond- 
ing rms noise fluctuations at high Galactic latitudes are in the 
range 60 ^ cr ^ 20 mK (Fig. [7]). We provide clean maps as 
the final product from our data processing but also maps with- 
out RFI elimination. To allow estimates on possible residual 
stray-radiation eff'ects it is also possible to extract the corrections 
which were applied. 
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